The inverse iron-bias in action in Abell 2028 
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Abstract 

Recent work based on a global measurement of the ICM properties find evidence for an increase of the 
iron abundance in galaxy clusters with temperature around 2-4 keV. We have undertaken a study of the 
metal distribution in nearby clusters in this temperature range, aiming at resolving spatially the metal 
content of the ICM. The XMM observation of the first object of the sample, the cluster Abell 2028, reveals 
a complex structure of the cluster over scale of ~ 300 kpc, showing an interaction between two sub-clusters 
in a "cometary" configuration. We show that a naive one-component fit for the core of Abell 2028 returns 
a biased high metallicity. This is due to the inverse iron-bias, which is not related to the presence in the 
spectrum of both Fe-L and Fe-K emission lines but to the behavior of the fitting code in shaping the Fe-L 
complex of a one temperature component to adjust to the multi-temperature structure of the projected 
spectrum. 
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1. Introduction 

The X-ray determination of elemental abundances in the 
hot gas of clusters of galaxies is in principle robust be- 
cause equivalent widths of the observed emission lines 
can be directly converted into the corresponding elemen- 
tal abundances (e.g., Mushotzky et al. 1996). Unfortu- 
nately the intrinsic simplicity of the measurement has 
faced the limitations of the X-ray satellites in terms of 
the shape of the instrumental response, bandpass, spec- 
tral and spatial resolution. The correct modeling of the 
temperature structure is crucial, in particular when deal- 
ing with spectra extracted from a large aperture centered 
on the core of clusters and groups, where photons coming 
from regions of different temperature and abundances 
are mixed together, given the presence of strong oppo- 
site gradients (i.e. cooler regions are more metal rich) in 
the temperature and metallicity profiles (e.g., De Grandi 
& Molendi 2001). 

The first important recognition of a bias in the measure- 
ment of elemental abundances was the description of the 
"iron bias" (Buote 2000). The previously found signif- 



icant subsolar values for the iron abundance in groups 
and elliptical galaxies stem from fitting a multi temper- 
ature plasma with a simple single temperature model, 
resulting in best-fitting elemental abundance biased low. 
Another trend in the ICM abundance versus cluster tem- 
perature has been shown with increasing evidence in re- 
cent years. Using the ASCA archive observations of 273 
objects Baumgartner et al. (2005) showed that clusters 
with gas temperature between 2 and 4 keV have a typ- 
ical mean abundance that is larger by up to a factor 
of 3 than hotter systems. Rasia et al. (2008) made a 
fundamental step in recognizing a possible bias analyz- 
ing mock XMM observations of simulated galaxy clusters 
and finding a systematic overestimate of iron for systems 
in the 2-3 keV range. They explained it as due to the 
fact that temperatures in the range of 3 keV are very 
close to the transition between the relative importance 
of the lines (Fe-L or Fe-K) used in determining the global 
iron content. Projection and low resolution effects can 
create a complex temperature structure, averaging in the 
same region different temperatures with different contri- 
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Fig. 1. EPIC image of A 2028 in the 0.5-2 keV energy band. The 
black circle indicate the core region used to extract the spectra of 
Fig. 2, discussed in the text. 
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butions of Fe-L and Fe-K emission. Simionescu et al. 
(2009) supported this explanation by analyzing a deep 
XMM observation of the high luminosity cluster Hydra A 
finding a biased high Fe abundance at a level of 35% in 
its core, dubbing for the first time this kind of bias the 
"inverse" Fe bias. 

With the aim of shedding more light on these issues by 
going beyond single aperture measurements, we started 
with the present work a study of the the metal dis- 
tribution in iron-richest systems of Baumgartner et al. 
(2005), selecting the 9 objects with Z > 0.6Zq (Anders 
& Grevesse 1989 units). Here we present the analysis of 
an XMM observation of one of these objects, Abell 2028. 

2. The XMM observation of Abell 2028 

The XMM observation of A 2028 clearly reveals two 
merging subclusters, a more luminous South- West com- 
ponent and a less luminous North-East one (see Fig.l). 
The "cometary" shape of the two subclusters seems 
to suggest the direction of the merging, with the SW 
main subcluster coming from the NW direction and the 
smaller subcluster coming from SE. The temperature 
map obtained with a single temperature model reveals 
that both subclusters have a cool core with a tempera- 
ture of 3 keV for the main subcluster and 2 keV for the 
smaller subcluster and they reach a peak temperature of 
~ 5 keV and ~ 4 keV respectively. 

If we focus on the spectrum of the core of the main sub- 
cluster, corresponding to the black circle in Fig.l, which 
is a circular region of 0.5' radius, corresponding to 44 
kpc. Even in this highest S/N spectrum the Fe-K line 
complex is not evident but still a single temperature fit 
returns a high metallicity (see Table 1). The IT model 
does not provide a good fit to the data, which clearly in- 
dicate the need for a range of temperatures from ~2 keV 
to 4-5 keV, and overestimates the abundance due to the 
inverse- Fe bias (see Tab.l and Fig. 2). The inverse- Fe bias 
is mainly driven by the behavior of the the Fe-L complex 



Fig. 2. MOS 1 spectrum of the core. With the solid green line the 
best fit model and residuals of the IT fit obtained in the broad 
0.5-10 keV band are shown; with the red line the fit with IT 
model in the 0.5-3 keV band; with the blue line the fit with the 
IT model with the parameters obtained by the 0.5-10 keV fit but 
with the abundance value kept fixed at the value obtained in the 
0.5-3 keV band. 



strength which falls rapidly with increasing temperature 
above roughly 3 keV: in a multi-T spectrum where the 
average T is near 3 keV, but there are temperature com- 
ponents below and above 3 keV, most of the the Fe L 
emission of this multi-T spectrum comes from the lower- 
T components in the spectrum. However, when one fits 
a (wrong) single-T model, the single temperature will be 
sufficiently high that the Fe L lines will be weaker, and 
thus the model will have to compensate by increasing 
the Fe abundance above the true value. 



Table 1. Results of the fits using IT and 2T models in different 
energy bands for the EPIC spectra of the core of A 2028. 



Model/Band 




Z 




xVdof 


(keV) 


(keV) 




(keV) 




IT 0.5-10 


3.0 ±0.2 


0.76 ±0.15 




252/168 


IT 0.5-3 


2.0 ±0.2 


0.33 ±0.06 




179/148 


2T 0.5-10 


5.0 (fixed) 


0.46 ±0.07 


1.4 ±0.1 


206/167 
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